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Abstract 

We solve the nonequilibrium dynamics of qubits or quantum spin 
chains (s=l/2) modeled by an anisotropic XY Hamiltonian, when the 
initial condition is prepared as a spatially inhomogeneous state of the 
magnetization. Infinite systems are studied analytically, yielding a uni- 
versal relaxation driven by quantum fluctuations. Particular long-lived 
excitations are found, for which the relaxation time diverges as a con- 
sequence of constructive quantum interference at degenerate stationary 
points. Those states are intrinsically entangled and may be of interest 
for performing quantum computation. We also numerically analyze finite 
samples to assess the extent of size effects. 

1 Introduction 

Magnetic quantum systems which displays a slow relaxation in the nonequilib- 
rium dynamics are good candidates for magnetic data storage devices. Low di- 
mensional magnetic systems, including magnetic chains and magnetic molecules 
in molecule-based magnets, may allow for a substantial reduction of the size of 
memory units to perform classical computation jT]. In addition, systems that 
develop quantum coherence in the nanometer scale may be of interest for the 
fabrication of quantum computers. Quantum computation is believed, at least 
in principle, to be much faster than classical computation, and its accomplish- 
ment rests in the ability of a quantum system to be in a superposition of many 
quantum states {entanglement). This property is the basis of the so called quan- 
tum parallelism, which presents a lot of advantages in relation to its classical 
counterpart, the one that is realized in a classical computer by connecting a 
large number of processors for parallel operation. Still, there is another impor- 
tant requirement for doing efficient quantum computation, a condition that is 
related to quantum programming or quantum software: we need robust Bell-like 
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states, or stable entangled states devised through constructive interference from 
the many quantum states available. The coherence time of those states must be 
much larger than the time employed in 1-qubit or 2-qubit operations 

The present contribution is related to this latter issue. It has been shown 
that quantum chains which are modeled through the Heisenberg Hamiltonian 
(or variants such as the quantum XY model), are intrinsically entangled The 
quantum XY model can be practically realized in low-capacitance Josephson 
junction arrays 0]; it also appears as an effective Hamiltonian for the interaction 
of quantum dot spins |S] , and it is also relevant for other condensed matter and 
QED systems which have been proposed as quantum computers [2]. One im- 
portant role of the XY model in quantum computation is that it can be used to 
construct a swap gate, which means that the time evolution can be consider as a 
source of entanglement 0. Moreover, an interesting connection between quan- 
tum phase transitions and entanglement has been recently proposed, suggesting 
that entanglement plays an important role in quantum critical phenomena [H]. 
Concerning the anisotropic XY model, a quantum phase transition is found as 
a function of the anisotropy parameter and transverse field. The latter may be 
envisioned as a quantum metal-insulator transition with the formation of a gap 
for the insulating phase. 

Here, we are able to construct analytic solutions of the anisotropic XY model 
in the thermodynamic limit, which display a very slow relaxation after being 
prepared in a non-homogeneous initial state. This is a peculiar nonequilibrium 
problem, since our quantum system is closed (does not exchange energy with 
the environment except when preparing the initial condition), and its relaxation 
to the homogeneous state is exclusively driven by quantum fluctuations and 
interference effects. Nonequilibrium properties of quantum systems have been 
recently investigated by several authors following different approaches [ZllHlinij 
but a complete theory is still lacking. To test our solutions, we also study 
finite systems through numerical calculations. Finite systems are interesting for 
practical applications, like the ones mentioned above. An important question 
to answer is whether a finite system will show similar relaxation properties, at 
least during finite periods of time, and how to monitor the magnitude of such 
periods. 

In spite that we solve a well defined model, some properties of the phe- 
nomenon seem to be universal for one dimensional systems. The key point 
addressed in this paper is how to prepare special long-lived superposition states 
in quantum magnetic chains, whose coherence is not wiped out by the strong 
destructive interference that appear at very long time. We call those states as 
critical, for reasons that will be developed in the next sections. 

A general nonhomogeneous initial state in a quantum spin chain can be 
Fourier decomposed. We then study the nonequilibrium dynamics through the 
time evolution of the Fourier component of the magnetization < Sg This 
approach has several advantages in spite that Sg is not an observable: one can 
probe the relaxation properties as a function of the wave vector Q, which is 
an important parameter in the theory, and for the particular case of the XY 
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model (and variants), the relaxation of the longitudinal component < Sq >t is 
independent of the initial state, except for a complex scale factor, i.e. 

<S'q >t^<S'Q >o nQ,t) , (1) 

where is a real universal function that depends on parameters of the Hamil- 
tonian. We are interested in solving the anisotropic XY magnetic chain with 
spin 1/2 in the presence of a transverse magnetic field 

N 

H = 2nuj {rS^ + [(1 + 7) S^S^+, + (1 - 7) ^J^+i] } , (2) 

where Tiojj is the exchange integral, r = {ojq/ujj) is the external transverse 
field in units of the exchange constant, < 7 < 1 is the anisotropy parameter, 
N is the total number of spins in the chain, and S" — {1/2) a" are spin 1/2 
operators. Special cases of this model are: 

i) the isotropic XY model with a transverse field (7 = 0); 

ii) the Ising model with a transverse field (7 = 1); 

iii) no magnetic field (r = 0). 

The chain is prepared in an arbitrary nonhomogeneous initial state and its 
magnetization dynamics is calculated in exact form under very general assump- 
tions. In the thermodynamic limit {N — > 00), and for long times, the system 
displays irreversibility in the form of a power law relaxation. This irreversibility 
can be ascribed to quantum interference effects which encompass the system 
at asymptotically long times. In analogy with classical optics, we handle this 
phenomenon using the stationary phase method jlO| . As a general rule, most 
of the components interfere destructively (shadow regions), and the relaxation 
is dominated by the contribution of stationary points, where the interference is 
constructive (illuminated zones). It may happen that some stationary points 
become degenerate in the parameter space. In optics, we get surfaces where 
the intensity of light diverges (caustics). In our problem, we get critical states 
whose relaxation time asymptotically diverges, yielding a very slow relaxation. 
Those critical modes do not occur in all cases, and we present a summary of 
our findings below. 

Models i) and ii) were discussed in Ref. JT]. Also Berim and Cabrera 
[T^ and Berim, Berim and Cabrera ^31 have solved similar problems for sev- 
eral versions of the 1-dimensional XY model, which include dimerization, bond 
alternation, and staggered magnetic field. All the models, including the one pre- 
sented here, are amenable of exact analytic treatment by fermionization meth- 
ods involving generalized Jordan- Wigner transformations I14| . Concerning 
relaxation properties, all the models studied can be grouped into two families, 
whether the energy spectrum has a gap or not jl3| . Only gapped models de- 
velop critical modes with anomalously slow relaxation properties. It remains 
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the question whether the above phenomenon is associated to a conservation 
law It is then interesting to study the behavior of Hamiltonian jSJ, since 
the anisotropy in the XY exchange and the transverse magnetic field lift the 
spin rotational symmetry and the time reversal symmetry, respectively. More- 
over, one can monitor the spectral properties (gap or not) with 7 and r. The 
present study confirms the general prescription stated in Ref. J3] for the XY 
family, i. e. anomalous long-time tails in the relaxation are features of gapped 
one-dimensional models. In this contribution, we will show that this anomalous 
slowing down of the relaxation is due to both, constructive quantum interfer- 
ence at degenerate or critical stationary points and nesting properties of the 
spectrum. 



2 The calculation 

Our calculation follows the general trends described in detail in Ref. |12[ I13L 
I15| . Hamiltonian ((SJ, with periodic boundary conditions, is written in terms of 
fermion variables (cj, cj) using the Jordan- Wigner transformation: 

= L,ic,+c])/2, S] = L,{c,-c])/2t, 5f = c]c, - 1/2, (3) 

with Lj being the sign string that adjusts the anticommutation relations for the 
(cj>c]) 

1=1 

For the fermion Hamiltonian (c— Hamiltonian), we solve the c-cyclic problem, 
where one neglects a correction term that comes from the boundary, as it is 
usually done when studying the thermodynamic limit (iV 00) jl4j. The 
c-cyclic Hamiltonian is diagonalized in two steps: 

i) we take the Fourier transform 

= V TV ^ "^^^ ^"'^ ' 

with k = 2rmT/N {-N/2 <m< N/2 - 1). At this stage, the Hamiltonian 
is not diagonal since the anisotropy couples the modes (fc, — fc); 

ii) complete diagonalization is achieved with a Bogoljubov transformation 

= Ukhk + ivkb'i^ , ril = Ukb\ - ivkb-k , (5) 
with real coefficients {uk,Vk) and Hamiltonian 

n^Y. ^^kvUk + c , (6) 

k 
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where C is a constant term and the quasi-particle dispersion relation is 
given by 

Afe = 2w J ^^(7 sin kf + {r + cos kf . (7) 

In the particle-hole representation, > and the ground state is the 
vacuum for the r] operators. We note that the spectrum given by Q may 
present a gap as a function of the anisotropy 7 and transverse field r. 
Several examples are available: 

(a) for 7=1 (Ising with transverse field), the spectrum presents a gap 
at A; = TT of magnitude 

A = 2|l-r| (8) 

in units of the exchange constant, which shows that there is a critical 
value Vc — 1 for the transverse field; 

(b) when r — (no magnetic field), we get a gap at A: = 7r/2 whose value 
is 

A = 27 , (9) 
the isotropic XY model being gapless; 

(c) with magnetic field and anisotropy, when r < 1 — 7-^, the gap is 
positioned at A: = arccos (^—j^^^ 1 with value 




1 -72 ■ 

Otherwise, its position and value are given by (jSjl. 



(10) 



Once the ground state and excitations have been determined, we proceed to 
calculate the dynamical properties of the magnetization. This is done using the 
identity 

{A),^TT[p{t)A]^Tv[pm{t)] , (11) 

where A{t) is an operator in the Heisenberg picture. This way, Hll(l relates the 
temporal evolution of the average at any time with the initial state p{0). We will 
probe the dynamics of the longitudinal (parallel to the field) Fourier component 
of the magnetization 

(S^q)^ = Tr [p{t)S^] = Tr [p{0)S^{t)] , (12) 

with 

N 

Sq^T. ^|exp(zQj) . 
i=i 

Exact closed forms for (|12l) are obtained when one assumes that the initial den- 
sity matrix is a functional of only one spatial component of the spins, a situation 
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that can be achieved in practice using strong non-homogeneous magnetic fields 
along a given direction, to prepare the initial state |12[I13| . We will denote those 
possibihties as p{0,S'^), with /i = x,y,z labehng the possible field directions. 
The calculation follows the general trends described in Ref. ^1^]. The details 
will be given elsewhere ^Hl- Note that the calculation of the dynamics of the 
transverse components Sq and Sq is much more involved, and it is not clear 
that one can obtain closed analytic formulae. For the above components we 
only obtained numeric results in small systems. 

For the longitudinal component, we summarize the main steps of the calcu- 
lation below: 

(a) transform Sq to fermion operators to get its time dependence. This step 
leads to: 

k— — 7r 

(13) 

where p ^ k - Q/2 , q = k + Q/2, Q.^{k,Q) = Ap ± A, , and (u,w) 
are the coefficients in the transformation ||SJ). Note that in H13II there are 
contributions from processes that do not conserve the number of particles; 

(h ) average with the initial density matrix. This means that we have to calcu- 
late the averages {nliiq)^ , ^ry.pTy^,^^ , (vl^^-q) ^ , and (rjlpi^l^^ in ^ 
over the initial state; 

(c) take the thermodynamic limit — > oo, i.e. introduce an infinite number 
of degrees of freedom. As a consequence, summations over the fc-space are 
replaced by integrals over the Brillouin zone. Using the symmetry of fc, 
one can reduce the integration to the interval [0,7r]; 

(d) find asymptotic analytic expressions for very long times, t oo. 

Note the order of the limiting processes of the two last steps. Indeed, for 
finite A^, the limit t ^ (x docs not exist. We will comment on this further in 
the following. For the intended scope of this paper, we will just consider the 
behavior of < Sq >t for the initial ensemble p{0, S^). All the cases, along with 
complete formulae, are given in |16| . We quote the exact result which comes 
from taking the continuous limit in H13|) : 

{Sq)^ = ^{S^q),£ dk {hik,Q) cos[n-{k,Q)t]+ (14) 
+g{k,Q) cos[n+ {k,Q) t]} , 



6 



where the frequencies 17^ are defined above and the functions (/i, g) are given 
by: 



h{k,Q) = 1 + [^^ + 2rcosA: cosQ/2 + cos2fc (cos2Q/2+ 

+7^2 gin2 Q/2) - sin^ k (sin^ Q/2 + 7^ cos^ Q/2)] , 
g(fc,Q)EE2-/i(fc,Q) . 



(15) 



The frequencies fi* which appear in and p4l) . correspond to particle- 
particle, hole-hole and particle-hole excitation processes that contribute to the 
time evolution of Sq. Note that for the regime t — > 00, the integral H14I) is 
strongly oscillatory, and integration over ordinary points leads to cancellations 
effects (destructive interference). The asymptotic behavior, obtained through 
the stationary phase method |l()j . is dominated by the contributions of sta- 
tionary points, where we get constructive interference. The detailed study of 
stationary points of fi* depends on 7, r and Q. The points fc = and k = ±7r 
are always solutions. In addition, the other stationary points are given as roots 
of a polynomial of fifth degree of the variable y — cosk. This means that the 
general problem has no analytic solutions, unless a fortuitous factorization is 
achieved. This way, we have only obtained analytic expressions for two impor- 
tant particular cases: model ii), the Ising model with transverse field (7 = 1 and 
arbitrary r); and model iii), no transverse field (r = 0) and arbitrary anisotropy 
7. As a general rule, the contribution of stationary points leads to an asymptotic 
relaxation with time in the form of a power law. We write 

n n 

where (n, ±) labels the stationary point of fl^ respectively. The number of 
stationary points and their degeneracy depend only on the parameters of the 
Hamiltonian (7, r) and the wave vector Q. The dominant time dependence of 
the (t) functions in l|16() is given by 

Jt{t)=Kt exp(ze±t) (J^^ " , (17) 

using the same notation as above. The important quantities in (|17|l are the 
relaxation rate and the exponent , which together determine the relative 
speed of the relaxation process. For the cases we calculated analytically, they 
are functions of (7,(5) or {r,Q). In Fig.l, we display the relaxation times as 
function of 7 for a given value of the wave vector Q (with r = 0). The asymptotic 
behavior is dominated by the contribution of five stationary points of $7+ : two of 
them are related to ti (the points fc = 0, tt), a third one (fc = 7r/2) is associated 
with T2, and from the two remaining we get the same relaxation time T3. The 
times T2 and T3 become asymptotically large in the vicinity of a critical value 
of 7, which is indicated in the figure, where the corresponding stationary points 
become degenerate. Some general comments are in order: 
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i) the asymptotic behavior oscillates with a small number of frequencies 8^, 
which are functions of (7, Q) or (r, Q); 

a) we get a whole set of exponents coming from the different stationary 
points, but the relaxation for t ^ co is driven by the smallest exponent 

in) depending on their variation in the parameter space, some stationary 
points may become degenerate. This fact drastically affects the relaxation. 
As remarked before, the r's asymptotically diverge in the neighborhood 
of critical points. Exactly at the critical values, one has to go some steps 
further in the asymptotic expansion (as many as the order of the degen- 
erate point), yielding a discontinuous change of the vo exponent. In the 
example of Fig. 1, the exponent i^q changes from 1/2 to 1/4 at the critical 
7(^ , signaling a slowing down of the relaxation process; 

iv) the above instance occurs at the loci of the critical curves Q = Qc in the 
parameter space, with Qc given by: 

, . / 2 arccosr , r < 1, , 

Qc{l = l.r) = ^ 2 arccos(i), r>l ' ^^^^ 

for the Ising case with transverse field, and 

Qc(7,?' = 0) = arccos f ^— ^ I , for < 7 < 1, (19) 
VI + 7/ 

for the anisotropic case and no field. Note that the gapless cases (7 = 
1, r = 1) and (7 = 0, r = 0) yield Qc —>■ 0, which is an ordinary stationary 
point. Also note that the pure Ising case (7 = l,r — 0) has no dynamics 
and has to be treated separately. 

The above results are exact in infinite chains, where the thermodynamic 
limit implies a continuous spectrum. In contrast, finite systems present quantum 
recurrences, with partial or total reconstruction of the initial state. For very 
small samples, we are able to identify a number of Rabi-like periods associated 
with the revival of the magnetization. Numerical computations show that those 
periods or quasi-periods become longer with increasing sizes, scaling almost 
linearly with A'^. This poses the question of how to define the time interval 
during which the finite system approximate the infinite system in a regular way, 
in the spirit of finite-size scaling methods. To estimate the relaxation, we have 
adopted the following procedure: we numerically analyze the dynamics of 
finite systems until the first partial reconstruction (this time depends on size 
and 7), and; b) we fit the time evolution to a power law similar to the one 
given by (|17l) . This is not an easy task, since there are several frequencies 
superimposed. A typical behavior of the magnetization is shown in Fig. 2, for 
A^ = 100 and for the same Q wave vector of Fig.l. We estimate the exponent 
from the envelopes of local maxima and minima, which are shown in the figure. 
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In Fig. 3, we display a summary of our findings for N — 100 spins. One observes 
a precursory slowing down of the relaxation around the critical value predicted 
for the infinite chain (see Fig. 1), with a typical rounding due to size effects. 
We comment those results further in the next section. 

3 Conclusions 

One interesting feature of the anisotropic XY model is the absence of any spin 
conservation law (we exclude the pure Ising limit without transverse field) . This 
fact makes its dynamics 'richer' than the one associated with other models, 
where the total spin or one of the total spin projections are conserved [T^IT^ . 
This appears as a source of entanglement, which is relevant for applications in 
quantum computation. Another peculiarity of the XY model is exhibited in ^ 
and H14|l . meaning that the dynamics of the longitudinal Fourier component of 
the magnetization is independent of the initial state, which enters just through 
a scale factor in the form of an average over the initial condition. This property 
is not shared by other models, like the Heisenberg XX Z model, where we found 
a strong dependence of the relaxation process on the initial state |17| . 

From our analytic and numeric computations, we were able to identify long- 
lived collective excitations, which have the character of dynamic spin density 
waves with incommensurate wave vector Qc- In the vicinity of Qc, we get 
t{Q) —^ oo asymptotically, thus reducing the damping of those excitations in 
relation to a general wave vector Q. Exactly at the point Q = Qc, the exponent 
Vq changes discontinuously (1/2 1/4) and the relaxation time r has a finite 
value. Note that critical modes only appear when the spectrum has a gap, in 
agreement with findings in other models of the XY family ^2 ^1 . 

Our theoretical results may be physically realized in a number of systems, 
from which the most promising to perform quantum computation are Joseph- 
son tunneling junction arrays. Quantum bits ('qubits') stored in low-capacitance 
Josephson junctions (JJ) simulate mesoscopic spins, and the coupling between 
JJ units mimic the spin-spin interactions that are usually found in several spin 
models 0]; in particular one can obtain an effective anisotropic XY coupling 
with a gap in the spectrum. Nonhomogeneous initial states considered in this 
paper, could be practically achieved in JJ systems by locally manipulating gate 
voltages, tunneling barriers or magnetic fields. This way, one may implement 
exact prescriptions to prepare those long-lived collective excitations as many 
times as necessary to perform specific tasks in a quantum computer JS] . Those 
states (which are not stationary states), are obtained by constructive interfer- 
ence from a macroscopically dense number of quantum states. This can be seen 
from the fact that the degeneracy of stationary points and the critical condition 
are associated to nesting of the one-particle spectrum around inflection points 
of the dispersion relation J^l- This peculiarity confers a degree of 'robustness' 
that is still present in finite systems, as shown in our example with — 100, 
and leads one to expect that quantum coherence will be maintained within 
operational levels when in contact with the 'outside world'. 
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FIGURE CAPTIONS 



Fig. 1 

The inverse of the relaxation times (adimensional units) associated with the 
stationary points of fl'^ as a function of 7 for the value Q = 27r/5. They were 
obtained through the asymptotic behavior at very long time {t — > 00) of the 
infinite chain. Here, the wave vector Q is fixed and the anisotropy is varied, 
crossing through the critical value = 0.528, which is indicated by an arrow. 
At 7(^, two branches are degenerate, and the corresponding relaxation times 
asymptotically diverge. Note that the branch of T3 does not exist for 7 > 7^- 

Fig.2 

Typical time evolution of the magnetization < Sq > (t) in a finite system, 
for a value of 7 away from criticality. The Q wave vector is the same of Fig.l, and 
the subindex n means that the magnetization is normalized to the initial value at 
t — 0. A regular relaxation behavior is obtained until ujjt « 250, after which we 
observed the onset of the first partial reconstruction. The exponent of relation 
ifTTjl is fitted using data from the upper (squares) and lower (circles) envelopes. 
In the displayed case, we obtained UQ{upper) « 0.476 and i^oilower) w 0.482, 
which are close to the predicted value vq = 1/2 for the infinite chain. 

Fig. 3 

The exponent 1^0 from numeric calculations for A = 100 spins as a function 
of the anisotropy parameter 7, for the same Q wave vector of the preceding 
figures. The analytic calculation for the infinite chain predicts a critical value 
7(^ = 0.528, where the exponent vq jumps from 1/2 to 1/4. The error bar 
represents the dispersion obtained from estimating the exponent using the upper 
and lower envelopes of the relaxation (see Fig. 2). 
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